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o. We show that AdS/QCD generates predictions for the rate of diffractive p-meson 

Oh electroproduction that are in agreement with data collected at the HERA electron- 



proton collider. 
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00 INTRODUCTION 

00 

o 

To date, the correspondence between string theory in five-dimensional anti-de Sitter 

O (AdS) space and four-dimensional QCD has enjoyed a number of successes (see [1-4] and ref- 

^ erences therein). In this letter, we demonstrate another success by showing that parameter- 

> 
\J~1 fr^^ AdS/QCD wavefunctions for the p meson lead to predictions for the rate of diffractive 

X 

H p meson production (7*^ — )■ pp) that agree with the data collected at the HERA ep collider. 

In previous papers [5, 6], we took a phenomenological approach and extracted the light- 
front wavefunctions of the p meson using the HERA data. We follow the same formalism 
here, except that we now use the AdS/QCD wavefunctions predicted in [7, 8]. 



THE ADS/QCD WAVEFUNCTION 

Brodsky and de Teramond have recently shown that, in what they call a first semiclassical 
approximation to light-front QCD [9], the meson wavefunction can be written in the following 



factorized form 



{x,Cv) = ^me''^^ (1) 



where L is the orbital quantum number and ( = ^/x(l — x)b {x is the hght-front longitudinal 
momentum fraction of the quark and b the quark-antiquark transverse separation). The 
function $(C) satisfies a Schrodinger-like wave equation 

d2 1-4L2 



d^2 4^2 +f^(C)J$(C) = M^$(C), (2) 

where ^(C) is the confining potential defined at equal light-front time. After identifying 
( with the co-ordinate in the fifth dimension, Eq. (2) describes the propagation of spin- J 
string modes, in which case U{() is determined by the choice for the dilaton field. We shall 
use the soft- wall model [10], in which 

U{0 = ^'C' + 2k\J-1). (3) 

This potential encodes the confinement dynamics of QCD and the challenge remains to 
derive it from first-principles QCD. 

Solving Eq. (2) with this potential results in eigenvalues 

M2 = 4fi;2(n + J/2 + L/2) , (4) 

which reproduces the correct meson mass spectrum. In particular, it predicts a massless 
pion {S = 0,n = 0,L = 0) and M^ = 2k^ for the p meson (5 = l,n = 0,L = 0). The 
corresponding eigenf unctions are [11] 

$(C) = «:x/2Cexp (-^\ • (5) 



It remains to specify the function f{x) in Eq. (1). This can be done by comparing the 
expressions for the pion EM form factor obtained in the light-front formalism and in AdS 
space [7] and it results in 

f{x) = M^x{l-x) . (6) 

The resulting wavefunction is thus 



^{x,0 = AT^^/^il^exp ( -'^^ , (7) 



where A/" is a normalisation constant. Assuming the meson is dominated by its leading qq 
Fock component, M is fixed by 

C^ dr 

d2bdx|0(x,C)P= / —^ -J\x) = l. (8) 

Jo x{l - x) 

Brodsky and de Teramond also have a prescription to account for non-zero quark masses 
[8]: A Fourier transform to /c-space gives 

1 / M2-\ 

0(a;, k) oc , exp --^ , (9) 

^x{l -x) V 2/t^ J 

where the invariant mass squared of the qq pair is M^^ = A;^/(x(l —a;)). For massive quarks, 
the invariant mass should rather be M^^ = (A;^ + mp(a;(l — x)). After substituting this into 
the wavefunction and Fourier transforming back to transverse position space, one obtains 
the final form of the AdS/QCD wavefunction: 

^(.,C) = iV^V;(r^exp(-^)exp(-^-,^) , (10) 

and A^ is fixed by the generalization of Eq. (8): 

r^^mexp(-^^)=K (11) 

Jo X[\ — X) \ K^X[\ — X) J 

This is rather similar to the original Boosted Gaussian (BG) wavefunction discussed in 

[12, 13] 



0^^(a;, C) oc x{l - x) exp ( ^^ ) exp ( - g^^/_^J exp [-^] ■ (12) 



If R^ = A/ H? then the two wavefunctions differ only by a factor of \/x{l — x), which is not 
surprising given that in both cases confinement is modelled by a harmonic oscillator. In 
what follows we shall consider a parameterization that accommodates both the AdS/QCD 
and the BG wavefunctions: 

0(x,C) oc [x{l-x)fe^^ ^-— j exp [- ^^,^^( _ ^^) ■ (13) 

COMPARING TO DATA, QCD SUM RULES AND THE LATTICE 

To compute the cross-section for 7*p — > pp we use the dipole model of high-energy 
scattering [14-17]. In this approach, the scattering amplitude is a convolution of the photon 



and vector meson qq wavefunctions with the total cross-section to scatter a qq dipole off a 
proton. QED is used to determine the photon wavefunction and the dipole cross-section can 
be extracted from the precise data on the deep-inelastic structure function F2. The details 
of this procedure can be found in [6, 13]. All that remains is to specify the wavefunction of 
the meson. 

The meson's light-front wavefunctions can be written in terms of the AdS/QCD wave- 
function 0(a;, C) [6]. For longitudinally polarized mesons: 

*t.(M) = ^V. (1 + ^1^) *.(..C) . (14) 

where V^ = \di, + dl and h (h) are the helicities of the quark (anti-quark). For transversely 
polarized mesons: 

where 6e*^ is the complex form of the transverse separation, b. Rather than using Eq. (11) 
to fix the normalization of we impose 

J2fd'hdx\^l,{b,x)\' = l, (16) 

h,h ^ 

where X = L,T. This means we allow for a polarization dependent normalization (hence 
the subscripts on (I)l,t)- For longitudinal polarization, the difference between the two nor- 
malization prescriptions leads only to sub-leading differences ~ M^/Q^ or m'l/Mp in the 
electroproduction scattering amplitude, where Q^ is the photon virtuality. For transverse 
polarization, the two prescriptions lead to a slightly different normalization but this can 
be attributed to the ignoring of higher Fock components in the wavefunction (since all of 
the normalization integrals are in any case only unity up to corrections due to higher Fock 
components). 

To compare with the HERA data we have to fix a value for the quark mass. We take 
rrif = 140 MeV, which is the value used in the fits to the deep-inelastic structure function, 
F2(x,g2) [6]. We use the CGC[0.74] dipole model [19, 20] (see [6] for details), ahhough the 
predicitions do not vary much if we use other models that fit the HERA F2 data [21, 22]. 
We also take k = Mp/y/2 = 0.55 GeV, which is the AdS/QCD prediction. 

Figures 1 and 2 compare the AdS/QCD predictions, shown as the solid blue curves, with 
the HERA data. The agreement is good and the disagreement at high Q^ is not unexpected. 



This is the region where perturbative evolution of the wavefunction will be relevant and the 
AdS/QCD wavefunction we use is clearly not able to describe that. It should be stressed 
that these AdS/QCD predictions are parameter-free. 

We are also able to compute the electronic decay width Fg+e-, which is related to the 
decay constant via 



ip 



47raL 



J I 
Using 



3r,+,-M,^ '/' 



we obtain Fe+e- = 6.66 keV, which is to be compared with the measured value Fe+e- = 
7.04 ± 0.06 keV [23]. 

Figure 3 shows the x^ per data point in the (/3, k) parameter space (see Eq. (13)) [28]. It 
confirms that the AdS/QCD prediction lies impressively close to the minimum in x^. The 
best fit has a x^ per data point equal to 114/76 and is achieved with k = 0.56 GeV and 
P = 0.47. Better fits to the data are possible, e.g. if one allows the parameters /3 and k, to 
depend on the polarization of the meson. However, given that we have not attempted to 
quantify the theory uncertainty we regard these as good fits. The curves resulting from the 
best fit are shown as the red dashed curves in Figures 1 and 2. 

We have previously shown that the twist-2 Distribution Amplitude (DA) can be related 
to (f)L{x, according to 

We note that J da: i^ix^ fi — ?■ oo) = 1 recovers the decay constant constraint. To compare to 

predictions using QCD Sum Rules [24] and from the lattice [25], we can also compute the 

moment: 

•1 

dx {2x- l)V(a;,/i) . (19) 



'0 

We obtain a value of 0.228 for the AdS/QCD wavefunction, which is to be compared with 
the Sum Rule result of 0.24 ± 0.02 at ;U = 3 GeV [24] and the lattice resuh of 0.24 ± 0.04 at 
/i = 2 GeV [25]. The AdS/QCD wavefunction neglects the perturbatively known evolution 
with the scale fi and should be viewed as a parametrization of the DA at some low scale 
/i ~ 1 GeV. Viewed as such, the agreement is good. 
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FIG. 1: Comparison the HERA cross-section data [26, 27]. Solid blue curve is the AdS/QCD 
prediction and the dashed red curve is the best fit. 
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FIG. 2: Comparison to the HERA data on the longitudinal to transverse cross-section ratio [26, 27]. 
Solid blue curve is the AdS/QCD prediction and the dashed red curve is the best fit. 
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FIG. 3: The x^ distribution in the (/3, k) parameter space. The AdS/QCD prediction is the white 
star. 
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